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We compute, at the first order in the fine structure constant, the parameters of 
the electromagnetic Lagrangian for the inhomogeneous Larkin-Ovchinnikov-Fulde- 
04 I Ferrell phase in Quantum Chromodynamics (QCD) and in condensed matter. In 



o 
o 



X 



particular we compute for QCD with two flavors the dielectric and the magnetic 
permeability tensors, and for condensed matter superconductors the penetration 



Qh! depth of external magnetic fields. 



I. INTRODUCTION 

The aim of this paper is to compute, at the first order in the fine structure constant, the 
parameters of the electromagnetic Lagrangian for the inhomogeneous Larkin-Ovchinnikov- 
Fulde-Ferrel (LOFF) phase in QCD and in condensed matter. In particular we compute 
in 2 Flavor QCD the dielectric constant and the magnetic permeability. For the condensed 
matter superconductor we compute the penetration depth of an external magnetic field 
at T = 0. The inhomogeneous superconductive LOFF phase was introduced in the con- 
text of ordinary superconductors forty years ago. In the original papers this phase 
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was discussed for weak ferromagnetic materials with an exchange interactions produced by 
the presence of paramagnetic impurities. It was shown that the spin splitting generates a 
separation of the Fermi surfaces. This separation, here denoted as 5fi, is proportional in 
metallic superconductors to the magnetic field H. For large enough 5fi, beyond the so called 
Clogston-Chandrasekhar limit ~ A /v^2 (A the gap for the homogeneous BCS case) 0, J], 
it can be energetically favorable for two electrons to form a pair with non vanishing total 
momentum |2q|. The main effect is an inhomogeneous gap with a space modulation which 
depends on its plane wave decomposition. The LOFF phase can exist also in QCD as a par- 
ticular realization of color superconductivity at non asymptotic densities, due to difference 
in Fermi momenta, as arising from different quark masses and from (3 equilibrium in dense 
quark matter. Translation and rotation invariance are broken and the space dependence of 
the order parameter may be that of a crystal. Such a crystalline phase of QCD might occur 
in compact stars, and suggestions exist that it may explain the variation patterns in the 
pulsars rotation periods (glitches) [^j]. 

The range of densities where the LOFF phase might be energetically favored is still matter 
of debate. In a recent paper [|| the intermediate density region has been studied and the 
possible spatially uniform candidate phases have been examined. The conclusion is in favor 
of a gapless CFL phase (gCFL) and, for immediately lower densities, in favor of the LOFF 
phase, based on the indications of the calculation in Ref. \[\. These results, if confirmed, 
would make more likely the occurrence of the LOFF phase at the pre-asymptotic densities 
of compact stars. 

An important point is the form of the condensate. Recent analyses (3, 0] point to cubic 
structures as the energetically favored form of the condensate. They are the body-centered 
cube (bcc) and the face-centered cube (fee), obtained summing 6 or 8 plane waves pointing to 
the faces or the vertices of a cube. The bcc structure seems the dominant one for 5 a near the 
CMo, limit 0; for ,ar g er va,ues of,, the fee structure b favored, flfl. Therefore here 
we consider three different structures, i.e. the one plane wave case (Fulde-Ferrell phase) and 
the two cubic structures. We perform our study in a well defined approximation, not based 
on the Ginzburg-Landau approach, but valid for A not too small |7|. This approximation 
is based on a convenient average over the sites of the crystalline structure defined by the 
condensate. The result can be described in terms of a multi-valued gap function possessing 
P branches, where P is the number of plane waves defining the crystal. Each of these 
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branches corresponds to a gap kA, k = 1, • • • , P with A the constant gap factor appearing 
in the LOFF condensate. 

In Sections |H] and IIHI we review the formalism employed to describe the LOFF phase 
and the results obtained by the approximation of ref . pj . In Section IIVI we discuss the 
LOFF phase in QCD. We discuss the problem of the Meissner mass and the determination 
of the parameters (dielectric constant and magnetic permeability) of the Lagrangian for the 
electromagnetic field. Differently from the homogeneous two flavor case (2SC) |9(, we find 
a correction not only for the dielectric constant, but also for the magnetic permeability. 
Since in QCD there is a rotated U{1) that is conserved Q|, this implies a constraint on our 
calculation scheme since the Meissner mass must vanish. We use this result in Section ITVl 
where we consider the LOFF phase in condensed matter. In this case the Meissner mass does 
not vanish and in general the magnetic field H should be expelled. This has been discussed 
in the LOFF superconductor Q within the Ginzburg-Landau approximation for the case of 
a gap with a space modulation ~ A cos 2qz and in |2j for the one-plane-wave case. In this 
Section we consider other crystalline structures and the region near the Clogston limit, far 
away from the second order transition point. 



II. GENERAL FORMALISM 



In this Section we briefly review the formalism we employ to describe inhomogeneous 
color superconductivity in QCD; modifications for condensed matter applications will be 
discussed in Section ITVl We consider QCD with two massless quarks having different chem- 
ical potentials /ii and fi2 arid we suppose that <5/i = — /x 2 1 /2 is slightly larger than the 
Clogston-Chandrasekhar limit Ao/a/2, where Ao is the value of the gap for the homo- 
geneous BCS phase. We work in zero temperature high quark density limit, which means 
that /i = (/ii + /i2)/2 3> Sfi. In these hypotheses the system can be supposed to be in the 
LOFF phase characterized by the following pattern of condensation: 

(Via C ^> = A J2 e 2 ~e Q/33 e fi (1) 

m=l 

where a,/3 are color indices, i,j are flavor indices and 2gn m is the total momentum of the 
Cooper pair. We will consider below three cases. The first is the one-plane wave Fulde 
Ferrel state with P = 1 (and n along the z-axis). In the second case we take P = 6 with the 
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six unit vectors n m pointing to the six faces of a cube (bcc). Finally we consider the case 
P = 8 with the eight unit vectors n m pointing to the eight vertices of a cube (fee). With 
the choice of phases for the plane waves as in Eq. the symmetry of the condensate both 
for P = 6 and P = 8 corresponds to the cube group. 

The reason to discuss only these cases is based on the results of p| . Here it is shown that 
the bcc is the energetically favored structure in the Sfi interval (0.707A — 0.95A ), while 
for 5/j, in the interval (0.95A — 1.12A ) the fee dominates. The apj 
is based on the so called High Density Effective Theory (HDET) 
and on an averaging procedure of the original Lagrangian over a region of the size of the 
lattice cell. In the HDET formulation one decomposes the fermion momentum p M in its hard 
part /if M and a residual momentum i^, i.e. p = fiv + £ where v M = (0, v); v is the fermion 
velocity and one neglects in £ the transverse part writing £ = (£ ,£\\v), with £\\ = £ ■ v. 
Since /ii ^ yU 2 we have two velocities here, but one can prove that in the large \x limit, 
vi = — v 2 + 0{b\ij\i). The momentum decomposition allows to define velocity-dependent 
fields, whose Fourier transform depends on £. The averaging procedure substitutes the 
inhomogeneous gap A(r) oc exp(2ign m • r) with a function of £ and v; the whole approach 
is justified if the velocity dependent fields are slowly varying over regions of the order of the 
attice size. Therefore this Lagrangian can only describe soft momenta. For more details see 
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Let us write the fermion propagator in this approach. Since we have four degrees of 
freedom (two flavors and two colors) we can use a compact notation introducing a base of 
velocity-dependent fermion fields ipA with A = 1, ... ,4. In this base the quark propagator 
assumes the form 

1 



with A,B = 1, 



Sab(v, 
, 4. Here V = 



1 V-15 AB 



D(v, 

;i,v),y 



Aab ^ 
y -A AB V-15 AB ) 

(1, — v), the gap matrix is 

(n n n i\ 



(2) 



A AB = A E (v,£ ) 



1 
0-10 
0-100 
10 



(3) 



and 



D(vJ) = V-£V-£-A 2 E (v,£ ). 



(4) 
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The effective gap is given by 



p 



As(v,4>)= E A e// (vn m ,£ ) , (5) 

with 

A e// (vn,e) = AO(E u )0(E d ) . (6) 

Here 

E u ,d = ±fy* =F gv ■ n + e (7) 

are the quasi-particle dispersion laws and e is the value of the energy at the pole of the 
propagator. 

III. 0(1) GAUGE INVARIANCE AND THE LOFF PHASE OF QCD 

The quark pair condensate ((H) breaks the electromagnetic group U(l) em since the pairs 
have total non-zero electric charge. There exists however a group which we call U(l) gen- 
erated by a linear combination of the electromagnetic charge Q and the T§(= As/2) color 
generator 

Q = Qf = Qij ® I aP - ® Tf (8) 

that remains unbroken as far as the transverse degrees of freedom are concerned. The 
residual symmetry embodied by Eq.fjHJ) implies a mixing angle 9 between the photon A and 
gluon G 8 : 

V3g 2 + e 2 

The in-medium vector potential fields A and G 8 are then given by 

Ai = -smeG^ + cosOAi , (10) 
= cos9 G 8 t+ sine Ai . (11) 

This phenomenon is similar to what happens in the two-flavor superconducting phase of 
QCD, the so-called 2SC model, where the quark pair condensate has the same color and 
flavor dependence of (0), but with q = 0. These results indeed do not depend on the space 
modulation of the condensate. Gauge invariance under U(l) implies that the polarization 
tensor Hij(p) of the Aj field vanishes for zero external momentum p = 0. To check this result 
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at the one loop level in the HDET we consider the self-energy and tadpole diagrams (see 



e.g. Fig.l in 
Lagrangians 



14j ) . They are computed by the propagator (J2J) and the effective interaction 



C 



qqA 



C 



qqAA 



V 



-v-AQ 



AB 



■v-AQ 



Xb 



AB 



( Q\b 
2fi+V-l 





2fi+v-e 



where e = ecos#, 



P^ = 9^-\(v it V v + V lt V v 



(12) 
(13) 

(14) 



and Q = 1/2 diag(+l, +1, -1, -1). 

For any crystalline structure we have the result 



iUijip) 



e 2 /i 2 

12 7T 3 

e 2 



dv r 2 V -£V- (£ + p) + V-£V- (l + p) -2A| 
AttJ iVj D(£)D(£ + p) 



— IdHp, 

67T 3 7 47T 



(a* + 'ii) s 



V 



+ (V - V) 



(15) 



(2/j + V -£)D(£) 

The former term on the r.h.s. is the contribution of the self-energy graph, the latter term is 
from the tadpole diagram. 

For the Fulde Ferrel one plane wave case we can distinguish two contributions, one from 
the pairing region, and the other from the blocking region. The pairing region is defined by 

Pi = {(^)|A = %(v.n )£ )| £= ^} . (16) 

The blocking region Vq corresponds to E u < or < or both. We have heuristically 
tested the vanishing of the Meissner mass in this case by computing numerically Ily(O) for 
the following values of the parameters /x = 400 MeV, 5/j, = 5/Ji = A /\/2 and taking the 
values of Sfi/q = 0.78 and of A = 0.24A that minimize the free energy Q. We find n(0) 
of the order of 10~ 3 , while the two individual contributions are of the order of 1, which 
points to an almost complete cancellation. This result was expected on the basis of gauge 
invariance and the rather accurate approximation employed to get (JZJ), i.e. dfi/fi <C 1. 

Next we consider the case of several plane waves with wave vectors with the same modulus 
q, but directed along the directions n m , m — 1, • ■ • , P. Here our formalism is based on the 
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approach discussed in Section |TTJ We have different regions Vk where pairing is possible. 
They are defined as follows 

V k = {(v, I A E (v, e) = kA} , k = 1, • • • , P (17) 

where Ae is in (JSJ) and one uses (J7|) with n — ► n m . This approximation is the result of an 
averaging procedure described in detail in [7j and is valid for £ of the order of q or smaller. 
Now in the self energy term in the r.h.s of (|TK|) the relevant contribution in the £ integration 
comes from the small £ region; therefore the approximation is adequate. On the other hand 
in the second term (tadpole contribution) the hard modes dominate and the approximation 
is no longer valid. Therefore we use U(l) gauge invariance to get information on the main 
features of the pairing regions for large £. Notice that the sum over k arises because plugging 
((HI) into (fTKJ) one has several terms, corresponding to different values assumed by the gap: 
kA (k = 1, • ■ • , P). Imposing gauge invariance one gets 

= v fl r d^dv k 2 A 2 
^\2Jv k 4tt (e + PA 2 ) 3 / 2 

+ f j dw ^ V* 2 + k ' 2A ' 2 ± g i (r , , r)) } an 

Jv k 47rVe + A; 2 A 2 ^ v^TfcW+(-2^ ^ i} ) J 1 J 

where is the region of the phase space where A^ = kA, but £ ~ //. Even if this equation 
does not determine Vk, it provides sufficient information to compute the Meissner mass 
in ordinary LOFF superconductors. An application of this result is discussed in the next 
Section. 

For small external momentum we get from Eq. (jl5j) and from the condition of vanishing 
Meissner mass, that 

iK(pk ~ -f^a \ A *iP 2 o + B ijP*P>\ ( 19 ) 



where we have defined 



^ f dHdv 2 , t Pf (Pidv 



+ 



4A| 



/J 2 (Z) /J 3 (/) 



(20) 



fc -> fc 47T 4 J /J 2 (/)' « fc f> fc 47T 

For the FF state we have two independent transverse tensors, IT^ = n n and Hj = H33 
while for both cubic structures we have 



n T (p) = \ L - ^ ) n«(p) (21) 
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as in the homogeneous case. The Lagrangian for the rotated photon can be written in the 
form 

where and Ajj can be obtained from Equations (f2T)|) . In general we have 

= (l + A o)t£^) ^ > V = [ l + ^> A o)Y^) • (23) 

The coefficients fj,gj assume different values according to the crystalline structure. At 
5fi = Sfii and Ao = 40 MeV we have, for the one-plane- wave (FF) fi = j'2 = +0.12, 
f 3 = +0.23, gx = g 2 = +0.31 x 10~ 2 and g 3 = +0.13 x 10~ 3 . For the body-centered- 
cube fbcc = +0.49 and gb cc = —0.09. For the face-centered-cube at dfM = 0.95A we have 
ffcc = +0.46 and gj cc = —0.09. For the 2SC case we have, in agreement with 0, Il7 |. 
/25c = 1 and g2sc = 0, showing absence of corrections for the magnetic permeability in the 
homogeneous phase. 

IV. PENETRATION DEPTH IN CONDENSED MATTER 

As an application of the result ()18jl . in this Section we give an estimate of the penetration 
depth of a weak static magnetic field in an ordinary condensed-matter LOFF superconductor 
at T = 0. We assume that the field is small enough to produce an exchange field and a 
sizeable paramagnetic effect so that the Fermi surfaces of the two pairing electrons are 
separated and the LOFF phase is formed (in particular the Clogston-Chandrasekhar limit is 
reached). At the same time we assume that the effect of the external field can be neglected 
in the gap equation. This is a strong assumption, as we know that the effect of an external 
magnetic field is a modulation of the gap, see e.g. 18]. As a matter of fact the paramagnetic 
effect, which is needed to produce the separation of the Fermi surfaces, and the diamagnetic 
effect, which is detrimental to superconductivity, are in general related. Therefore the 
original proposal of is now considered only as an ideal case. The actual experimental 
activity points to layered superconductors where one can minimize diamagnetic effects by 
choosing the external magnetic fie ld p arallel to the layer. We refer the interested reader to 
the specialized literature (see e.g. |19] and references therein) for a discussion. For the time 
being we study the effect of the magnetic field in the idealized case where its effects on the 
gap equations can be neglected. 
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Let us assume that a plane surface (yz plane) divides the space into two parts, one 
containing the superconductor in the LOFF phase (half-space with x > 0) and the other 
containing matter in the normal phase. At the interface the magnetic field H is parallel to 
the yz plane. We take H along the z— axis and the vector potential A directed along the 
y— axis; A depends only on x and we assume V • A = 0. 

In condensed matter the LOFF condensate has the form similar to ((TJ: 

p 

(?pCtp) = A(r) = A e 2iqn ™ r (24) 

m=l 

where the ip are non relativistic, two components spinor fields describing electrons. The 
Lagrangian can be written as follows [19I ] 

/ 



(V -l + ev- A) 5 ab + 
-A E (v,£ ) 5 ab 



5/2 a ab 



-A E (v,£ )5 ab 
(V ■ £ — ev • A) 5 ab + 5fia : 



ah 
3 



Xb 



(25) 



where we have used the effective Ag(v, £q) approximation; 8[i is proportional to the exchange 
field acting on the electron spin and the term Sfia^ describes a paramagnetic coupling. In 
f!25j) a, b are spin indices and \ a are Nambu-Gorkov fields. The Lagrangian includes the 
coupling to the external vector potential field A. The penetration depth is defined by j^j 



6 H Ju 



dx H(x) 



where H is the value of the magnetic field outside the superconductor. If rotationa 
metry holds (the BCS and the cubic structures) one gets, using previous hypotheses 

1 



(26) 



sym- 



5 



2 

7T ^0 



dp ■ 



p* - n T (p) 

where Ht(p) is computed by (j2"Tj) in the static po = approximation; for Ely (p) we have 

dv 



(27) 



n«(p 



where we have defined 



e 2 pfttiv 2 



4tt 2 



E 



J(C,(3,kA) = — 



e~PZ + k 2 A 2 



Vj J(C,A kA) + dUij , 



e + PZ + k 2 A' 



(28) 



(29) 



'(e- (3) 2 + k 2 A 2 ^(^ + (3) 2 + k 2 A 2 _ 

with (3 = p • v/2 and v is the direction of Fermi velocity v. The sum over k in ()28j) goes 
from k = to k = P. For the BCS homogeneous case there is one term: k — 1. For the 
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one-plane wave case there are two terms, one with k = corresponding to the blocking 
region, where A = 0, and the other one with k — 1, corresponding to the pairing region. 
For all the other cases the sum runs from k — 1 to k — P. As a matter of fact, within our 
approxlmatl0 n 3, for t ne s_ with mOT e t h a n one p.ane wave, pairing is possible in 
the whole phase space. In these cases, as discussed in Section IV. A of Q|, one can identify 
the blocking region with the domain where only the branch with gap PA of the dispersion 
law contributes. The two terms correspond to the two terms on the r.h.s. of Eq. ()15j) . 

We are interested in Type II superconductors where the relevant momenta in Eq. (|2*7)l 
are p ~ 0. For the BCS case one gets 

2 r°° 1 1 

5 w - / dp 2 = (30) 

7i Jo p z + m M m M 

where m 2 M = v is the squared Meissner mass in the BCS phase. 

Let us consider now the LOFF phase. For the one-plane-wave case we have 

IMP) = -m 2 M (*« - 3 l Q ^ e j+ lf d{ £ g v t 0, J (fc El* a) ) . (31) 

The dependence on the total momentum of the Cooper pair 2q is in the definition of V\. 
It is convenient to consider the tensor IT with q along the x-axis. The relation between the 
two tensors is (sum over k, I) 

ny(p) = R lk (6)R 3l (6)fl kl (p), (32) 

where Rij{0) is the rotation matrix which brings q along the x axis, that is along p (the 
direction of the gradient of the magnetic field). II has two independent components n n = 
II22 and II33 and the superconductor is characterized by two independent penetration lengths 
which we compute in the London limit. We do that for the following values of the parameters 
[l]: S/i = 5fii, z q = 0.78 and A = 0.24 A . In this case we get: 5i ~ 2.6 5 L and <5 3 ~ 1.4 5 L 
where 5l is the London penetration depth in the BCS case. Within our approximation, 
consisting in neglecting terms of the order 5fi/ fi our results are compatible with those of j^; 
for example for 63/61 we find agreement with [2] within 10%. 

We notice that for A — ► 0, near the second order phase transition, the pairing region 
vanishes, whereas the blocking region Vo is the whole Fermi surface. From ()31)1 we see that 
rijj(0) vanishes, and from ()27|) we get that both 6i diverge, which means that the FF is no 
longer a superconductor, in agreement with the result of jl(. 
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We can repeat the analysis for cubic crystalline structures (bcc and fee). In these cases 
one can exploit the residual discrete symmetry and only one penetration length is present. 

The transverse component of the polarization tensor is obtained by (J2HJ) using the appro- 
priate set of plane waves. In the London limit only p ~ are relevant and one gets: 

n T (0) = 2m M g4 ^? + &A* { W + M + + { ^ -°J 



* r djdw fc 2 A 2 

2^ L Anr IC2 1 U2 A 2"\3/2 " 



-m 



4vr (£ 2 + PA 2 ) 3 / 2 ' 



In the second line we have used the result expressed by Eq. ()18j). Numerically we get for 
the London penetration length in the bcc case 5 ~ 0.695^ at = <5//i; for the fee case we 
get 5 ~ 0.525l at 5/i = 0.95A , where, according to [3] there is a transition from the bcc to 
the fee LOFF phase. 



V. CONCLUSIONS 



We have used the high density effective theory formalism to compute the low energy prop- 
erties of the electromagnetic Lagrangian of the LOFF phase in QCD and condensed matter. 
We have shown that in QCD the rotated photon associated to the unbroken U(l) group is 
screened both electrically and magnetically. We have computed near the Chandrasekhar- 
Clogston point the dielectric tensor and the magnetic permeability tensor for the one-plane- 
wave, the body-centered cube and the face-centered cube crystalline structures. In condensed 
matter we have computed the penetration depth of an external magnetic field. In the Lon- 
don limit the penetration depth is proportional to the London penetration depth of the BCS 
case with coefficients that assume different values according to the crystalline structure. 
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